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Abstract 

We study the spectral zeta functions associated to the radial Schrodinger problem with 
potential V{x) = x^^'^ + ax^^~^ + (A^ — l/4)/a;^. Using the quantum Wronskian equation, 
we provide results such as closed-form evaluations for some of the second zeta functions 
i.e. the sum over the inverse eigenvalues squared. Also we discuss how our results can be 
used to derive relationships and identities involving special functions, using a particular 
5F4 hypergeometric series as an example. Our work is then extended to a class of related 
(— 1 1 PT-symmetric eigenvalue problems. Using the fused quantum Wronskian we give a simple 

0^ method for calculating the related spectral zeta functions. This method has a number of 

applications including the use of the ODE/IM correspondence to compute the (vacuum) 
nonlocal integrals of motion Gn which appear in an associated integrable quantum field 
theory. 



*E-mail: J.Watkins@kent.ac.uk 



1 



1 Introduction 



A common set of models to be studied in quantum mechanics are the anharmonic oscillators 
defined by the Schrodinger equation 



-tp" + x^^'4j = Ei;. (1.1) 

Any eigenfunction must obey the requirement that ■?/'—)■ as a; —;■ +oo and additionally satisfy 
either the Dirichlet or Neumann conditions at the origin. The anharmonic oscillators include 
the two exactly-solvable cases of the harmonic oscillator (M = 1) and infinite square well 
(M = oo). For these two models the eigenvalues can be found in closed form, whereas for 
general M numerical methods must be used. The Dirichlet eigenvalues {E^} and Neumann 
eigenvalues {E^} are associated to the spectral zeta function 

oo _ 

WKB estimates show that E"^ oc fc^z+i as A; — )■ oo and hence ^zp(s) converges for 3?(s) > 
(M + 1)/2M. In the two exactly-solvable cases Z^{s) can be given in terms of the Riemann 
zeta function. However Z^{s) can also be calculated for s G N, the most concise example being 



given by^ 



where M > 1 and 



0rr(l - (t(1 T 2)) 



a = . (1.4) 

M + 1 ^ ^ 

The zeta functions Z^(2) can be also be calculated [35], although these functions involve a 5F4 
hypergeometric series and no example of a closed form expression has been given. 

There has been much study related to these particular spectral zeta functions [35, 33, 34, 
16, 31], with one significant result being the determination of 'sum rules' which relate together 
different Zz^{s) when s G N [35]. An elegant example of this is found for the quartic anharmonic 
oscillator, the first few sum rules being 

Z+(l) =2Z„(1) 
2Z+(2) =Z_(2) + 3Z_(1)2 
2Z+(3) =9Z„(1)3 - Z„(l)2 - Z„(3). 

The techniques for obtaining the sum rules, which will be covered later, are the consequence 
of a wider study known as exact quantisation, a review of many such problems being found in 
[32]. For a review of the physical applications of spectral zeta functions in general, we direct 
the reader toward [22]. 

We will consider a generalisation of (1.1) given by the Schrodinger equation 



- < + \^x^^' + ax""'-' + ^—^j tP = Eij (1.6) 

with boundary conditions on the half-line. The eigenfunction criteria is that — )■ as a; — )■ +00, 
with the wavefunction having either of the behaviours tp^ ~ X2~^^ or tp^ ~ X2~^ as x — 0. 
Such eigenvalue problems have a long and rich history and have been investigated in many 



^Calculations for two related zeta functions were originally given in [35] and in [16]. 
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contexts such as quasi-exact solvability [30, 14], integrable models [21, 23, 9] and 'PT-symmetry 
[12, 13, 18, 17]. 

The two boundary behaviours ip- and ip^ respectively define two sets of spectra {E'^} and 
{E"^}, which we say are the eigenvalues of the 'regular' and 'irregular' radial problem[27]. To 
these spectra we associate the zeta functions Z^{s)^ which for A = 1/2 are the functions defined 
in (1.2). However for general A the regular and irregular problem are related by the analytic 
continuation A — )■ —A [21], giving the zeta-function identity 

Z+(s,A) = Z„(s,-A). (1.7) 

Throughout A will be restricted from taking the values 

A = ±^((2mi + l)(Af+l) + a) and A = ± (^ma + ^(M + 1)) (1.8) 

where mi,m2,m3 G Z"*", although when a = the condition becomes G Z"^ [21]. The 
first restriction excludes the possibility that any could be equal to zero and the second 
restriction is to ensure linear independence of the ip- and ips^ [21], a condition which is necessary 

2M 

for later work. Finally, as the eigenvalues are known to have large-fc behaviour E^ oc /c^^+i, 
the associated spectral zeta functions Z^{s) will converge when 3?(s) > (M + 1)/2M. As we 
will be studying Z^{s) for s > 1, we also restrict M > 1 throughout. 

This paper will be divided into two sections and will investigate the radial eigenvalue prob- 
lems (1.6). In Section 2 we use the well-known method involving Green's functions to calculate 
Z^:{s) when s G N. Then we will use the quantum Wronskian equation to derive appropriate 
sum rules and give examples where Z^{2) can be written in closed form, including examples 
involving the anharmonic oscillators. Similar techniques are then used to derive functional 
relations for two different hypergeometric series which appear in the expressions for ^zp(l) and 
Zip (2). These functional relations can recover many specific, known properties of hypergeomet- 
ric series, as well as providing new identities. In Section 3 we will extend our work to a class of 
eigenvalue problems commonly investigated in PT-symmetric quantum mechanics. The zeta 
functions Zx{s) of these problems are studied for s G N and found to be constructible by sum 
rules involving Z^{s). These expressions for Zk{s) are then used in conjunction with the ODE 
/ IM correspondence [8, 21, 19] to calculate the vacuum nonlocal integrals of motion G„ which 
appear in a related quantum integrable field theory [4, 5]. 

2 Radial Schrodinger problems 

Given that the eigenvalues {E^^} have no known closed form for general M, naively we might 
expect that associated spectral zeta functions Z^{s) cannot be given explicitly, except in the 
exactly-solvable cases. However when s G N the zeta functions are computed using the process 
given in [35]. Although Z^{s) can be defined by analytic continuation when s G Z~ [33], 
generally there is no known method to calculate their explicit forms when s takes general 
values. Therefore we elect to consider 2'=p(s) only when s G N, which we notate Z=p(n). 

For a Hermitian eigenvalue problem the eigenfunctions are complete and the Green's func- 
tion is written as 

niE;.y)^±til^, (2.1) 

fc=0 

where {E^} are the associated eigenvalues. As the irregular problem is generally non- Hermitian, 
the eigenfunctions are not guaranteed to be complete and consequently the Green's function 
cannot be written as in (2.1). Instead Z^{n) is defined by analytic continuation as in (1.7), 
which is valid except for A as in (1.8). 
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As the eigenfunctions in the Hermitian case are necessarily orthonormal, the zeta functions 
are calculated by 

Z_{n)= I R{0;xi,X2)R{0;x2,Xs} . . . R{0;Xn,Xi)dx, 

where M"*" denotes the integration over all positive space in n dimensions. To evaluate this 
repeated integral, R{0]x,x') is decomposed as a combination of ipL and tpR, two linearly- 
independent wavefunctions which solve the Schrodinger equation when E = 0. For a radial 
problem the two wavefunctions have specific asymptotic behaviours; ipi must obey the boundary 
conditions at the origin and we require that ipR — )■ as x — ?■ +00. The Green's function is 
written as 

R{0]x,x') = ^4jl{x<)4jr{x>) 

where x< = min(x,a;'), a;> = max{x,x') and W = ipRip'i — '4^'r4^Li the Wronskian of the two 
solutions. Thus 

^-H = ^ I ^L{xi)^R{xn) JJ %l)L{xi)^R{xi+i)dxidx2 ■ ■ ■ dx^- (2.2) 

Now we specialise to the regular radial problem (1.6), which is solved for = by 

<T-1 tT — 1 1^ 

tpL = X M_s^^a-x{'2(JX'') and 'ipR = x W-s^-ax{2ax''). 



Here M and W are the Whittaker functions [25] and the Wronskian of the two solutions is 
given by 

After making the change of variables p = 2ax{ and q = 2ax2 , we find that 

2-20-^/1 _|_ £0 _|_ roo 

Z4l) = J^,L'l^ / P'''~'W.^,.,Ap)M.^Mp)dp (2.3a) 



4-r(l + 2aA) JO 
and 

Z-m = 16^P(\ + La) / ?"-'W^%,-A.(g) I p'-'MK_.^,Mdpdq. (2.3b) 

The evaluation of these integrals requires known identities from (7.625) in [25]. The integral 
(2.3a) can be calculated immediately and found to converge for a > 1/2, reproducing the 
restriction M > 1 which was imposed earlier. Completing the integration in (2.3b) is more 
complicated and can be handled in a similar way to the calculation for the anharmonic operators 
as given in [35]. Completing the integrations gives the first two zeta functions 

cr^-^-r(i + ^ + aA)r(2a(l + A))r(2a) / 1 + ^ + gA, 2^(1 + A), 2a\ 



4-r(l + 2aA)r(i + ^ + cr(2 + A)) ' Vl + 2aA, i + ^ + a(2 + A) 

and 

^4-4<x - r(i + ^ + Aa + m)r(i + ^ + Aa + n) 1 



16<7 V{l + 2a\ + m)V{l + 2a\ + n)m\n\ im + n + I + 2ail + 

Lm,n=0 



xG' 



22/ 2 "'^^2 

33 



i-crA, i + ffA, l + m + n + ^ + 2a{2 + X) 



/ + m + n + cr(4 + A) - i, / + m + + a(4 + 3A) - |, 

(2.4b) 
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where {a)b is the Pochhammer symbol, pFg is the generahsed hypergeometric series and G is 
the Meijer-G function [25]. The latter two functions are both evaluated at x = 1. 

By setting a = 0, the Whittaker functions in (2.3a) and (2.3b) are replaced by Bessel 
functions and the integration techniques directly follow those in [35]. The zeta functions then 
take on the neater forms 

7 n m ^^'^'^r(a(i + A))r(a)r(i-a) 

" = °) = 4v/?r(l - a(l - A)) ('-'"^ 

and 

0F(t3-4^ r(2a(l + A))r(a(2 + A))r(2a) 



Z_{2,a = 0) 



4i+-^(l + A) r2(l + aA)r(l + a(2 + A)) 

/i + .A,2.(l + A),.(2 + A)2.,.(l + A) 



Setting a = in (2.4a) also recovers (2.5a) directly by Dixon's theorem. 
2.1 Radial sum rules and simplifications {a = 0) 

For M > 1 the growth rate of the eigenvalues is sufficient to ensure the convergence of the 
spectral determinants 

l-1^V (2.6) 



D_(E,M,a,A) = Z}_(0)n( 



Here -D-(O) is a regularising prefactor designed to vanish when some = 0. By analytic 
continuation we define D^{E, A) = D^{E, — A) except for the values of A given in (1.8). 

In addition to the analytic continuation relating D^{E) and D+{E), there is the quantum 
Wronskian equation 

2Xu^ = u^D_(uJE,-ia)D+{uE,ia)-u^D_{ujE,ia)D+{uJE,-ia), (2.7) 

where 

u = exp^iira). (2-8) 

The sets of eigenvalues {E^} and {E^}, having no relationship except at (1.8), are thus related 
to each other by a simple expression involving their associated spectral determinants. 

Two clarifications need to be made regarding the quantum Wronskian equation. First is 
to note that when a ^ the spectral determinants in (2.7) refer to differential equations 
with imaginary coupling constants on the term, a problem which is remedied in [32]. 

Secondly, for general a, the quantum Wronskian links together four different spectral determi- 
nants, whereas setting a = reduces this number to two. For this latter case the subsequent 
results are much less complicated and so a = is fixed for the remainder of this section. 

Using (2.7), we can derive sum rules relating together different Z^{n). For small E the 
spectral determinants can be written in terms of their zeta functions [33] as 

~ ZAn) 



D^{E) = D^{0) exp g -^E-j . (2.9) 
Using the technique of [35], (2.9) is substituted into (2.7) and the coefficients of the powers of 
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E are compared to obtain the radial sum rules. The first few of these are 

= A^iZ_(l) + iV_iZ+(l) (2.10a) 
= + iV_2^+(2) + {N^ - N2) (Z+(l) - Z_(l))2 (2.10b) 

= + iV_3^+(3) - ^ (iVs - N^m) (Z+(2) - Z_(2))(Z+(1) - Z_(l)) 

_ 1 (iV3 - + 2iVf) (Z+(l) - Z_(l))=^ (2.10c) 

where 

Na = sin(7rcr(A + a)) csc(7rcrA). 

Under certain parameter choices, (2.10) allows for Z^{n) to be given in terms of spectral zeta 
functions with a simpler form. As an example, the choice cr(A + 2) = m G N implies that 
= 0. At such a point (2.10b) does not include Z^(2) and hence, for m G N, we find 

+ ^ ' /V ^ 24-4'^r2(l-3(T + m)r2(l + (T-m)sin(47ra) ^ ^ 

This expression is much less complicated than the general form for Z^(2) given in (2.5b), 
reducing the 5F4 hypergeometric series to a product of gamma and trigonometric functions. 

Example 2.1. The anharmonic cubic oscillator has the zeta values 

8(V5- l)7r4 



Z+(2) 



5¥r^(|)P(|) 



and 



5¥r5(|)P(A) 5Sr8(|)r3(|) 5Sr(|) % ' 



In addition to finding closed form expressions for Z^{n), there is a reason to investigate the 
composite zeta functions 

Z{s) = Z+{s) + Z_{s) and Z{s) = Z+{s) - Z_{s). (2.13) 

For the potential V{x) = with eigenfunction condition ip G L2(]R) and eigenvalues {Ek}, 

Z{s) and Z{s) are respectively the 'full' and 'skew' zeta functions, defined as 

Z(s) = y-f- and Z{s) = y^-^. 

k=0 ^ k=0 ^ 

For the anharmonic oscillators, these definitions correspond to (2.13) as the spectrum {Ej.} is 
the interlacing union of {E^} and {E'^}. For general a and A this cannot be guaranteed. 

Although generally used as a convenient notation, simplifications in the explicit forms of 
Z{n) and Z{n) will be of use later. The radial sum rules are used to find reduced forms for 
Z(2) and Z{2) by realising points where N2 = ±A^_2- For m G Z+ we find 

,,,, 7rV^'^r2(l - 2cr)sec2(7ra)sec(27ra) , 



4i-2.r4(l - (7)r2(| - a - m)r2(i - a 



m) 



2 ^ .'"y- V2 

and 

7r5sec2(^)tan(^) 
Z 2,(T = l/4 = , ^ ^' ^ \\ . 2.15 

^ ' ' ' 64r4(|)r2(3±A)r2(3^) ^ > 

There are no simplifications for ^(2, A = 1/2) and no general simplifications for Z{2), as A 
would conflict with (1.8). 
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Example 2.2. The anharmonic sextic oscillator has the zeta value 



Z{2) 



32r4(|)r2(|)r2(|)' 



2.2 Functional relations between hypergeometric series {a = 0) 

While the radial sum rules in (2.10) are useful for determining properties of Z=p(n), they can 
also be used to determine properties of the special functions appearing in the explicit forms for 
Z^[n), as in (2.5). The are two principle reasons that such derivations can be made. First is 
the analytic continuation (1.7) and the second is that the radial sum rules can be written as 

Z+(1) = -^Z41) (2.16a) 

These identities implicitly specify functional relations for the special functions appearing in 
(2.5). By substituting (2.5a into (2.16a) and using (1.7), we recover the identity 

r(cT(l-A)) r(a(l + A)) sin(7rcT(l + A)) 



r(l - cr(l + A)) r(l - a(l - A)) sin(7ra(l - A)) ' 

which is effectively two copies of the Euler reflection formula for the gamma function. 

The next application is to the zeta functions Z=p(2), the explicit forms involving a 5-F4 
hypergeometric series as in (2.5b). Therefore we expect to derive properties of the function 

TiW^ w( i + «^A,2a(l + A),a(2 + A),2a,a(l + A) \ 
^ ^ - " ' Vl + ^A, 1 + 2aA, i + a(2 + A), 1 + a(l + A);* ' 

Substituting Z^{1) and Z^(2) into (2.16b) gives the functional relation 

sin(7rcr(A + 2))A-''^T{2a{l + A))r(a(2 + A)) 



sin(7ra(A - 2))(1 + A)r2(l + aX)Ta + a{2 + A)) 



-F(A) = (2.18) 



4-^r(2a(l-A))r(a(2-A)) ^, 4^^-'anlT\l - 2a)T\a{l + X)) 



(A-l)r2(l-CTA)r(| + cr(2-A)) ' ' r4(l-a)r2(l-CT(l-A))r(2(T)sin2(7ra; 

'sin^(7ro-(l + A)) sin(7ro-(2 + A)) 2 sin(7rcr(l + A)) sin(7rcrA) 



X 



sin2(7ra(l - A)) sin(7ra(2 - A)) sin(7ra(l - A)) sin(7ra(2 - A)) 



Furthermore, numerical testing indicates that this identity appears to hold when < 3ft((j) < 
3/4, although our result only applies when o" G M. 

By choosing a and A to take special values - precisely those which give simple forms for 
Z+(2) - we recover further properties of J-". For example the simplification (2.11) gives, for 
m G N, the identity 

, , (m-3a)r^((T)r2(l + 2(T-m)r2(l-2(T)r(i + 4(T-m)sin3(27r(T)csc2(37r(T) 

^ (2 — TTl (J) — 

^ ' ' 0F4™+i-4'^r2(l -3cT + m)r2(l + a-m)r(2CT)r(6(7-2m)r(4a-m)sin(47ra)' 

When m = 1 this identity is verified by Dixon's theorem but we have been unable to find 
any known identities accounting for general m G N. While T can be rewritten in terms of 
4F3 hypergeometric series [35], this appears to provide no extra information and the functional 
identity (2.18) is not apparent. 

The exact form of 2^- (3) is known to contain Appell series [35], functional relations for 
which being obtained by the above methods. We would then expect relationships between two 
different Appell series in terms of 5F4 hypergeometric series. 
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2.3 Functional relations between hypergeometric series {a ^ 0) 

The techniques used in Section 2.2 also apply when a 7^ 0. In this situation the sum rules found 
from (2.7) are more complicated, featuring four separate zeta functions. Substituting (2.9) into 
(2.7), the coefficient of E in the subsequent expansion gives, after the transformation a — ?■ ia, 

D_(0, a)£'+(0, -a) (u;^~^Z+(l, -a) + uj'^^^^^ Z_{\, a)) 

= D_{0, -a)D+{0, a) (cj^~^Z+(l, a) + cj^+^Z_(l, -a)) . (2.19) 



where -D=p(0) is given in [18] as 



D_(0, a) oc (2a)^-^-^/^ n^^^^^'^K. - (2-20) 

The explicit form for Z-{1) features a 3F2 hypergeometric series as in (2.4a). Therefore (2.19) 
is implicitly a four-term functional relation involving this function. Splitting into real and 
imaginary parts, (2.19) is actually two independent functional relations satisfied by the 3F2 
hypergeometric series. We define the function 

X) ^ I ^ /^l + f + 2a(l + A), 2.\ 

y[a, Aj - ^ 2a + f + aA)r(i - f - aXf^' [l + 2aX, l + ^ + a{2 + X))' 

where the hypergeometric series converges for 3fJ((T) < 1/2. Substituting (2.4a) and (2.20) into 
(2.19) gives, after taking real and imaginary parts, the pair of four-term relations 

^/ XX ^/ XX sin(7r(7(l-A))r(l + 2(7A)r(2(T(l-A)) , 
and 

C(r.\\ C(r.\\ cos(^^(l - A))r(l + 2aA)r(2a(l - A)) 

^(a. A) - G{-a, A) = ,,,(,,(1 + ;,))r(l - 2aA)r(2a(l + A)) "^^ " "^^^ " 

(2.22b) 

This pair of equations allows for the derivation of three-term equations featuring G(a, A). One 
such example is given by 

r(i + 2(7A)r(2(7(i-A))r(i-2(7(i + A)) , , , 

G{a, A) = ^ ' ^ ^ ^ (sin(27rfr)e;(«, -A) - sin(27raA)e?(-a, -A)) 

vri \L — la A ) 

(2.23) 

Now Qia.^ A) has been shown to satisfy a particular three term functional relation. From (2.23), 
specific parameter choices allow for identities which only feature ^(«, A) twice. One example is 
given by A = 1 — l/2cr. When this is the case Q{a^ A) simplifies to a hypergeometric series, 
which is in-turn expressed in terms of gamma functions by Gauss' theorem. The subsequent 
relation is that 



Y\\ 2a)r(2 2a) - e;(2(5, 1 - l/2a) + e;(-2(5, 1 - 1/2(7). (2.24) 



r(2a)r(l - a{\ + 5))r(l - o{\ - 5))r(2 - 4a) 



Similar such reductions can be found by any simplification of the 3F2 hypergeometric series in 
(2.21) being substituted into the three-term relation in (2.23). 
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Figure 3.1: Stokes rays for M = 2.7 and K = 1,2,3 are given by the solid (blue) lines. The 
dashed (red) lines are the relevant anti-Stokes rays and we choose the quantisation contour 
to asymptotically tend towards these, while avoiding the branch cut along the imaginary axis. 
The different boundary conditions thus define separate spectra. 




3 PT-symmetric Schrodinger problems 

Now we consider a class of non-Hermitian eigenvalue problems related to the radial problems 
of the previous section, given for i^' G N by the Schrodinger equation 

- 0" + (^i-lfitx)''' - a{ixY'-' + = E<p. (3.1) 

The boundary conditions are that (j) G L^(C(x)) where C{x) is a quantisation contour that 
asymptotes towards the anti-Stokes rays with complex arguments — 7r/2 ± tt{K + 1)/ (2M + 2). 
These anti-Stokes rays are shown in Figure 3.1. Such a contour must also avoid the branch 
cut along the positive imaginary axis, employed to ensure that the potential is single-valued 
for general M . 

Surprisingly these Schrodinger problems, with complex potentials and boundary conditions 
defined in the complex plane, can have an entirely real spectrum [12, 13, 18, 28]. The suggested 
reason for the possibility of spectral reality lies in the PT-symmetry of the problem, meaning 
that (3.1) and the boundary conditions are invariant under a reflection in the imaginary axis. 
Specifically the spectrum is known to be real if the "PT-symmetry is unbroken, meaning that 
the eigenfunctions are invariant under the combined VT operation [13]. The revelation of a 
non-Hermitian problem exhibiting real eigenvalues has since spawned a huge research effort 
with hundreds of papers. Recent review articles of the subject are also available [10, 19]. 

Rather than investigating spectral reality, our interest will be in the relationships between 
the radial problems in the previous section and the PT-symmetric of this section. For fixed 
M, different choices of K will result in quantisation contours which asymptotically lie within 
different pairs of Stokes sectors. Generally these different problems are then expected to have 
unique spectra for fixed M, a,A, as was demonstrated in [13]. To these separate spectra we 
associate the eigenvalues {Ej^} and the spectral functions 

CKiE)^CKiO)ll(^-^) ZKis)^Y.7^s- (3-2) 

k=0 ^ ^ fc=0 ^ ^' ^ 

Despite the different spectra appearing to have no obvious analytic relationship, in [29, 21] 
a number of results are given which link together different Ck{E) and hence the associated 
eigenvalues. Following work in [21], using the normalisations in [17], one such result is the 
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'fused' quantum Wronskian equation 

- a;(^'+i)^Z}„(u;^+^E,2^+i«)Z}+(cZ;^+iE, (3.3) 

where M > 1 and A is again restricted to exclude the values in (1.8). 

The fused quantum Wronskian is essentially a simple extension of (2.7), which can be 
verified by setting K = and Cq{E) = 1 in (3.3). The function Cq{E) is understood to 
be a constant as the choice K = Q dictates that the quantisation contour of the eigenvalue 
problem asymptotically lies within two adjacent Stokes sectors, meaning an empty spectrum 
[21]. Therefore any sum rules derived from (3.3) must reproduce those derived from (2.7) by 
setting K = Q and defining Zo{s) = 0. 

In contrast to (2.7), the fused quantum Wronskian will feature only two spectral determi- 
nants on the right hand side when a 7^ and K is odd. At such values of K the sum rules 
must be of similar form to those in (2.10). 

3.1 Fused sum rules and simplifications 

Calculating Zxin) directly can be handled by the Green's function method used previously. 
This was first implemented in [26, 15] for the anharmonic (a = 0, A = 1/2) oscillators, where the 
results were used to verify conjectures on the spectral reality of the "PT-symmetric problems. 
However there is an easier method than direct calculation, using sum rules and the exact forms 
for Zz^{n) calculated in (2.4) and (2.5). To determine sum rules between Zxin) and Z=p(n), 
Ck{—E) is expanded near the origin as 

Ck{-E) = Ck{0) exp ^^(-Er^ . (3.4) 

The fused sum rules are calculated by substituting (2.9) and (3.4) into (3.3) and comparing 
coefficients of the powers of E. Temporarily suppressing the dependence on K when a ^ 0, 
the first few fused sum rules are 

Zk{1) = -L1Z4I) - L_iZ+(l) (3.5a) 
Zk{2) = L2Z42) + L_2^+(2) + {Lj - L2) Z{lf (3.5b) 

Zk{3) = -L:,Z43) - L^sZ^S) + ^ (^3 - L2L,) Z{2)Z{1) 

+ ^{L,-3L2L, + 2Ll)Z{lf (3.5c) 

where 

La = sm{iTa{K + 1)(A + a)) csc{7ra{K + 1)A). 

For a = the fused sum rules apply for all K, whereas for a ^ this is only true when K is 
odd. In the latter situation Zx{'n,a) is given in terms of Z^{n, —a) for K = An — 3 and in 
terms of Zzp{n, a) for K = 4n — 1. 

To calculate Zxin) it is now only necessary to calculate Z_(ra), . . . , Z_(l), determine 
Zj^{n, A) from (1.7) and then substitute these into the appropriate fused sum rule. This method 
implies that Zxin) is always real- valued, despite no guarantee that the eigenvalues will be real. 
The explanation is that all PT-symmetric systems have a real characteristic equation, implying 
that the eigenvalues are either real or appear in complex conjugate pairs [11]. 

Given the fundamentally differences between the radial problems and the "PT-symmetric 
problems, there is a priori no reason to expect such simple connections between the differing sets 
of spectra. Specific examples of the fused sum rules can be written in many ways, depending 
on the relationships between Z_{n) and Z_^.{n) as determined by the radial sum rules (2.10). 
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Example 3.1. The anharmonic quartic oscillator has the sum rules 

= 2Z„(1) 
Zi(2) = Z_(l)^-Z_(2) 
2Zi(3) = 3Z_(2)Z_(1) + Z4lf. 

As with the radial sum rules, the fused sum rules are also useful for deriving simplifications of 
Zx{n). For example, comparing (2.10b) and (3.5b) shows for m G N that 



4^4-4(Tr2 



T\l-2a) csc2(7r(T) 



Z^(2, a = m/(A + 2), a = 0) = ^ ^ _ ^^y.^, _ ^ ^^y.^, _ x 

(csc(37rcr) + csc(7r(j))^ sin(47ro-(A' + 1))(1 + sec(27rcr)) 



csc^{TTa{K + 1)) sm\27ra{K + 1)) (1 + 2 cos(27r(T))2 sin^(7r(T) sm{2TTa{K + 1)) 
Example 3.2. The anharmonic cubic oscillator has the zeta value 



(3.7) 



V5j 5ir4(|)n(|)- 

Certain choices of M and K reduce the quantum Wronskian to the form Ck{—E) oc D^{±E)D^{±E) 
and hence give the zeta function identity Zxin) = (+l)"Z(n). Additionally (3.5) will allow 
for the form Zxin) = Z{n) + f{Z^{n — 1), ... , Z=p(l)) under appropriate parameter choices. 
Therefore the simplifications for Z{2) and Z{2), as in (2.14) and (2.15), apply to Zk{2). Two 
examples of these are : 

Example 3.3. The quartic oscillator has the zeta value 

Z2(2,« = 0,A = 3/2)= f^Vr^ 



.2/ V3, 

Example 3.4. The anharmonic sextic oscillator has the zeta value 

(3 - 2V2)tx^ 



Z2{2) 



I6r4(f)r2(|)p 



3.2 Nonlocal integrals of motion 

An application for the fused sum rules in (3.5) is within the 'ODE/IM Correspondence' [20, 8, 
21,9] which links the eigenvalue problems in this paper to a class of integrable models. The 
correspondence, which was established by functional relations such as the quantum Wronskian, 
is summarised in-depth in the review article [19]. 

Two essential components for realising the correspondence are the continuum analogues of 
the T- and Q-operators which appear in the study of integrable systems such as the six-vertex 
model. These analogues were introduced in [4, 5] to study the integrable structure of conformal 
field theory. When M > 1 and a = the vacuum eigenvalues of these analogues, denoted T(s) 
and Q{s), are related simply to the spectral determinants of the radial and PT-symmetric 
problems. This remarkable connection was shown in [8, 21] and is summarised by the relations 

T{s) = Ci{-us^) (3.8) 

and 

Q{s) = ^f-i'^^") (3-9) 
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under the parameter identifications 



P 



aX 



and u 



a 



2a-2 



a- 



(3.10) 



Considering (3.9), direct relations are given in [5] between the nonlocal integrals of motion 
(IMs) Hn, used in a power series expansion for Q{s), and the zeta functions Z^{n). Additionally 
conjectured in [5] were analytic continuations for Z_((2n — l)/{2a — 2)) in terms of local IMs 
/2n-i and for Z_[—n/a) in terms of nonlocal IMs Hn- Therefore there appear to be good 
reasons to consider the spectral zeta functions in relation to the nonlocal IMs. Specifically we 
will work with Gn, which were given in [4] to be nonlocal IMs appearing in the power series 
expansion 



T(s) = T(0) + ^G„s2". 



(3.11) 



n=l 



These functions are found directly from the integral 



Gr, 



2n 



dui 



Ul 



n 

j>i <- 



4 sin 



dvi 



sm 



ru2 

/ dV2. ■ 




Jo 




Jo 








)] 


n 


2 sin 










with the first of these being given simply as 

47r2rfl 

Gi 



2/32 



r(i -/32-2p)r(i-/32 + 2p)' 



(3.12) 



(3.13) 



The nested integration required to calculate Gn is reminiscent of the nested integration required 
to compute Z_{n) as in Section 2. The difference is that while Z_{2) can be written exactly, 
no general form for G2 is known, although alternative approaches have been used in [24, 6, 7]. 
To calculate Gn, (3.8) is expanded on both sides as a power series by using (3.4) and (3.11). 
Comparing powers of s and using T(0) = Ci(0) = 2 cos(7rcrA) [4, 21], the first few G„ are found: 



T\l - a) cos(7raA)Zi(l) 



/n-\ 2(7-2 

G,^2(|) - 

/rr\ 4o-— 4 

G2 = ( - ) r^(l - a) cosiivaX) {Z,{lf - Z,{2)) 



2 

6 (2 



G3=l f^V" ' r^(l - (t) cos(7raA) {Zi{lf + 2Zi{3) - 3Zi(2)Zi(l)) 



(3.14a) 
(3.14b) 
(3.14c) 



Now that Gn are expressed in terms of the zeta functions Zi{n), they can be written in terms 
of Z^{n) by (2.10). This circumvents the need to calculate G„ directly, expressing them instead 
in terms of zeta functions which can be calculated from work in the previous sections. Given 
that Z-{n) becomes more complicated as n increases, identities for Gn are expected to reflect 
this property. However the integral expression for Z-{n) in (2.2) can always be given in terms 
of some infinite power series. This provides a computational advantage over (3.12), which has 
not been directly evaluated for n > 2. 

The vahdity of (3.14) can be checked by first substituting (2.5a) into (3.5a). The general 
identity for Zi{l) is then substituted into (3.14a) and, after the variable changes in (3.10), found 
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to agree with (3.13) exactly. The process is simply extended to G„ and the higher nonlocal 
IMs giving, as an example, that 

47r^r2(l -2cr)sec(7raA) cos^(7ra) 



r2(l - (t(1 - A))r2(l - cr(l + A)) V amHmil - A)) smHn<7(l + A)) 

\2) ^ \ sin(27r(7A) ' sin(27r(TA) ^'j ^ ' 

where the identifications in (3.10) must be taken into account. 

This expression for G2 can be checked to agree with some special cases given in [6]. Practi- 
cally (3.15) provides little extra information over the methods given in [24, 6, 7] due to the lack 
of known analytic properties of the 5F4 hypergeometric terms appearing in (2.5b). However the 
simplifications for Zi^(2), for example (3.7), allow for concise special values of G2 to be found, 
in some cases being as simple as the general formula for G\. 



Example 3.5. For = 2/5 and p = 11/10 the second nonlocal IM is given by 



45r4(|) 

The above methods can be extended to compute G3 and higher nonlocal IMs. Simplifications 
available by the sum rules would mean that, for specific parameter choices, G3 can be given 
without having to calculate a general form for Z^{3) directly as in Section 2. 



4 Summary and future work 

We have calculated the spectral zeta functions Z^{1) and Z^{2) pertaining to the radial prob- 
lems in (1.6). Using the quantum Wronskian, sum rules were established which give relation- 
ships between the different Z;p(n). As well as illustrating the deep connections between the 
two problems, these sum rules allow for the simplification of Z^{2) with many examples found 
in the anharmonic oscillators. The physical reasons for such simplifications, if any, is not yet 
understood. 

Our methods were then applied to determine special properties of the gamma function 
and specific hypergeometric series. We have not found another way to obtain the functional 
relations between these hypergeometric series and a separate derivation would be of interest. 
This technique could be applied in general to any Schrodinger problems where sum rules can 
be established between the zeta functions. Furthermore the functional relations established 
encodes many known, closed-form identities for specific hypergeometric series. A full study of 
this topic may therefore prove interesting in the future. 

In Section 3 the zeta functions Zxin) of a related, PT-symmetric problem were calculated 
not directly but using the fused sum rules relating them to Z^{n). This requires fewer cal- 
culations than the direct method and serves to highlight the deep connections between the 
real-line and monodromy problems. After calculating Z^ln), we showed how the ODE / IM 
correspondence can be used to calculate the vacuum nonlocal IMs Gn- This was restricted to 
a = 0, although we expect our work will be useful in testing the correspondence postulated in 
[9, 1, 2], where this restriction is not imposed. 

We deliberately excluded the possibility of a zero-energy eigenvalue. However our initial 
investigations show that this need not be the case. In particular we have found that spectral 
zeta functions can be used to detect the presence (and classify the order) of zero-energy excep- 
tional points in non-Hermitian eigenvalue problems. A natural application of this method is to 
the PT-symmetric problems specified in (3.1), where the fused sum rules are used. We hope 
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to work more on this in the future. 
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